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ABSTRACT 



We show that the anomalous retrograde perihelion precession of Saturn Am, 
recently estimated by different teams of astronomers by processing ranging data 
from the Cassini spacecraft and amounting to some milliarcseconds per century, 
can be explained in terms of a localized, distant body X, not yet directly discov- 
ered. From the determination of its tidal parameter K, = GMx/r x as a function 
of its ecliptic longitude Ax and latitude fix, we calculate the distance at which X 
may exist for different values of its mass, ranging from the size of Mars to that 
of the Sun. The minimum distance would occur for X located perpendicularly 
to the ecliptic, while the maximum distance is for X lying in the ecliptic. We 
find for rock-ice planets of the size of Mars and the Earth that they would be 
at about 80-150 au, respectively, while a Jupiter-sized gaseous giant would be at 
approximately 1 kau. A typical brown dwarf would be located at about 4 kau, 
while an object with the mass of the Sun would be at approximately 10 kau, so 
that it could not be Nemesis for which a solar mass and a heliocentric distance 
of about 88 kau are predicted. If X was directed towards a specific direction, 
i.e. that of the Galactic Center, it would mimick the action of a recently pro- 
posed form of the External Field Effect (EFE) in the framework of the MOdified 
Newtonian Dynamics (MOND). 

Subject headings: Experimental studies of gravity Experimental tests of gravitational 
theories 
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PACS: 04.80.-y,04.80.Cc 



1. Introduction 



Anderson fc Nietol (120091 ) recently examined some still unexplained a nomalies connected 
with astrometric data in the solar system. They a re the flyby anomaly (jAnderson fc Nieto 
20081 ) , the Pioneer anomaly (jAnderson et al.l Il998l ) , the secular change of the Astronomical 
Unit (IKrasinskv fe Brumberd 120041 ) and the increase in the eccentricity of the Moon's orbit 
( [Williams fc Boggsll2008l ). In fact, there is the possibility tha t also a fifth anomaly does 
actually exist: the anomalous perihelion precession of Saturn (Ilorioll2009ah . 



The corrections Am to the standard Newtonian/Einsteinian secular precession of the 
l ongitude of t he perihelion^ m of Saturn, estima ted with the latest versions of the EPM 
( Pitjeva 2009 ) and INPOP (jFienga et al. 2009a ) ephemerides by including some years of 
continuous radiometri c ranging data to Cassini in additi on to data of several types spanning 
the last century, ar<| (lPitievall2008l ; iFienga et ai1l2009ri ) 



Ari7p it = —6 ± 2 mas cy , 
ArijFie = — 10 ± 8 mas cy" 1 ; 



(1) 
(2) 



both are non-zero at a statistically significant level (3<r and 1.2a, respectively) and they are 
compatible each other since their difference is equal to 4 ± 10 mas cy" 1 . At the momenlR 
no corr ections Am estimated with the DE ephemerides by NASA JPL are available. Ilorio 
(j2009al ) unsuccessfully examined several possible dynamical explanations in terms of both 
mundane, standard Newtonian/relativistic gravitational physics and of modified models 
of gravity. Anyway, further analyses of extended data sets from Cassini with different 
dynamical force models are required to firmly establish the existence of the anomalous 
perihelion precession of Saturn as a genuine physical effect. 

Here we will show that the existence of a localized distant body (planet X/Nemesis), 
modeled in neither EPM nor INPOP ephemerides, is a good candidate to explain a secular 
perihelion precession of Saturn having the characteristics of eq. ([TJ-eq. (T5]): indeed, 



1 m = u + Q, where u is the argum ent of perihelion and Q is the longitude of the ascending 
node, is a "dogleg" angle (jRoyll2005l ). 

2 The formal, statistical error in the Pitjeva's result is 0.7 mas cy -1 ; Pitjeva (E.V. Pitjeva, 
private communication, 2008) warns that the realistic uncertainty may be up to 10 times 
larger. 

3 Plans to impr ove the ephemerides of Saturn with VLBA observations to Cassini exist 
(I Jones et al.ll2009h . 
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contrary to a massive ring usually adopted to model the action of the minor asteroids 
and of the Trans Neptunian Objects (TNOs), it yields a retrograde secular perihelion 
precessions and the constraints on its distance for different postulated values of its mass are 
consistent with several theo retical predictions put fo rth to accommodate some features of 
the Edgeworth-Kuiper belt (ILykawka fc Mukaill2008l ). Concerning Nemesis, it would be an 
undiscovered stellar companion of the Sun which, moving along a highly elliptical orbitfl 
would periodically disturb the Oort cloud being responsible o f the periodicity of about 2 6 
Myr in extinction rates on the Earth over the last 250 Myr (IWhitmire fc Jackson! Il984l ; 



Davis et al.lll984l ); the Nemesis hypothesis has also been us ed to explain the measu remen ts 



of the ages of 155 lunar spherules from the Apollo 14 site (iMullerl |2002| ) . See e.g. Ilorio 
( I2009bl ) for further details. Interestingly, such a proposed explanation of the anomalous 
perihelion precession of Saturn in terms of pointlike dark matter is, to a certain extent, 
to be considere d as degenerate since also the MOdified Newtonian Dynamics (MOND) 



( IMilgromlll983l ) predicts certain subtle effects in the planetary region of the solar system 



that may mimic the a ction of a distant mas^f) located in the direction of the Galactic Center 
(GC) (lMilgromll200ll ). 



2. The action of a distant body and of MOND's External Field Effect on the 

perihelion of a planet 

A hypothetical, still undiscovered body X, located at heliocentric distance rx>r along 
a direction hx, where r is the distance of a generic known planet P of the solar system, 
would impart on it a perturbin g acceleratio n Ax consisting of an "elastic" Hooke-like term 
plus a term directed along fix (jloridl2009bl ) 



where 



A x ~ — /Cr + 3/C (f • hx) hx, 
r ._GMx 

= 5 



(3) 



(4) 



is the so-called tidal parameter of X. Note that th e acceleration of eq. ([3]) derives from the 
following quadrupolar potential flHogg et al.lll99ll ) 



K 



r 2 — 3(f • hx) 2 ] 



(5) 



4 It should have semimajor axis a = 88 kau and eccentricity e > 0.9 (IWhitmire fc Jackson 



1984). 



5 It should not be confused with the supermassive black hole in Sgr A* whose action on 
the solar system's planets is, as we will see, at present undetectable. 
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Ioriol (j2009bl ) worked out the orbital effects of eq. on the longitude of perihelion w 
of a planet P by means of the standard Gauss perturbing approach with the assumption 
that fx can be considered constant during an orbital revolution of P. More specifically, 
the Gauss equation for the variation of the longitude of perihelion w of a planet under 
the action a small p erturbing acceleration A, whatever its physical origin may be, 
flBertotti et alJl2003h 



is 



ZD 



V 

nae 



-A r cos / + A T ( 1 + - ) sin / 



2 sin 2 



5 in, 



(6) 



where A r , A T are the radial and transverse components of A, respectively 
semimajor axis, eccentricity and inclination, respectively 



a, e 



I are the 

^GM/c 3 



T) = \A - e 2 , n = ^uivi/a" is 
the unperturbed Keplerian mean motion, p = a(l — e 2 ) is the semi-latus rectum, and / is 
the true anomaly counted positive anticlockwis e from the perihelion. The Gauss variation 
equation for the node fl is flBertotti et al.l 120031 ) 



n 



na sin 



-.A, 



sin(cj + /), 



(7) 



where A v is the normal component of A. It must be recalled that, in order to make 
meaningful comparisons with the estimated corrections Am to the standard perihelion 
rates, they have been obtained by processing the planetary data in the standard ICRF 
frame which is a frame with the origin in the (known) solar system's barycenter and having 
the mean ecli ptic at J2000 ep och as reference plane with the x axis directed towards the 
Vernal point (I Johnston! Il999l ) . This is particularly important when there is some physical 
feature, like a static body in a given direction as in our case, which breaks the spatial 
symmetry: assuming that fix, which is actually a-priori unknown, coincides with one of 
the frame's axes just to simplify the calculation is, in principle, incorrect. The analytical 
calculations for a non-privileged direction of fix are very cumbersome^; at the end, one is 
left with an expression of the kind 



where F is of the form 



(w) x = JCF(a, e, /, Q, u; A x , Px), 



F = J2G t (a,e)T t (I,Q,uj;Xx,f3x) 



(9) 



in which Gi are complicated functions of the semimajor axis a and the eccentricity e of P, 
while Tj are trigonometric functions of the inclination /, the longitude of the ascending node 
Q and the argument of perihelion u of the planet P perturbed by X, whose ecliptic longitude 
and latitude are Ax and fix- Releasing such analytic expressions would be, actually, 



6 The software MATHEMATICA has been used. 
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extremely space-consuming and of little help: they have to be numerically computed for 
given values of a, e,I, As a result, by comparing eq. (jSJ) to the estimated correction 
Am for a given planet P like Saturn, one has the tidal parameter K, of X as a function of 
Ax, Px- 

At this point, it is interesting to note that the form assumed by the Extern al Field 
Effec#l (EFE) in the planetary regions of the solar system in the recent study by iMilgroml 
( 120091 ) has exactly the same functional dependence of eq. ([3]), provided that 



K 



A 9 

n 



n — > r t 



GM & 



6.833 kau, 



(10) 



where ( IBegeman and Broeilslll99ll ) A = 1.27 x 10 10 m s 2 is the characteristic acceleration 
scale of MOND, and 

n x -> -x. (11) 

Indeed, ( Milgrom 2009 ) uses a frame with one coordinate axis directed toward^ GC; as we 
will see below, it is approximately the opposite of the x axis of ICRI0 Although such an 
effect would manifest itself in the strong-field regime existing in the planetary regions of the 
solar system, the functional form of q depends on the form of the MONDian interpolating 
functional /i(X) in the transition region in which X ~ 1, i.e. approximately at r t . 

Let us start to examine just the case of a dark object placed in the same direction of 
GC. The right ascension a and declination 5 of GC, assumed coincident with Sgr A*, are 
flReid fc Brunthalerlliooi ) 



a GC = 17 h 45 m 40.045 s , 
S GC = -29°0'28.l". 

The relations among a and 5 an d the eclip tical longitude A and latitude f3 are, from 
standard spherical trigonometry (jRoyll2005l ). 



(12) 
(13) 



sin j3 
cos A cos (3 
sin A cos (3 



cos e sin 5 — sin a cos 5 sin e 

cos a cos 5 = —0.054, 

sin e sin 5 + sin a cos 5 cos e 



-0.097, 



-0.99, 



(14) 
(15) 
(16) 



7 It is one of the non-linear features of MOND according to which the internal dynamics 
of small system s, like the solar system, does depend on the external field of a larger system 
S, like the Milky Way, in which s is embedded. 

8 The Galactic external field is, indeed, the source of the centripetal acceleration A c m A 
of the Sun during its motion of revolution around GC. 

9 It is the z axis in the frame used by IMilgroml (120091 ) . 

10 X is the ratio of the total gravitational acceleration felt by a body to A . 
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where e = 23.43 deg is the obliquity of the Earth's equator to the ecliptic. Thus, since, by 
definition, 0° < A < 360° and -90° < (3 < +90°, eq. fl33}-eq. flE) yield for GC 



Agc 

Pgc 



183.15°, 
-5.61° 



(17) 
(18) 



which tell us that GC is approximately directed in the opposite direction of the ICRF x 
axis. The GC longitude and latitude yield for the tidal parameter JC of a hypothetical 
X/Nemesis object 



/C Pit = (2.1 ±0.6) x 10" 26 s" 
/C Fic = (3.5 ±2.8) x 1(T 26 s" 

or, equivalently, for the MOND quadrupole parameter —q 

0.34 ±0.10, 



<?Pit 



0.6 ±0.4. 



(19) 
(20) 

(21) 
(22) 



Note that iMilgroml (120091 ) predicts 10 2 < —q < 0.3 for the relevant range of values for the 
Galactic field at the Sun's location, and for a variety of interpolating functions; thus, eq. 
fl2U)-eq. (|22D tell us that 

-g>0.2; (23) 

the upper bound is less tight being 



0.4 - 1. 



(24) 



It must be noted that IMilgroml (120091 ) made certain simplifications in his calculations that 
should be taken into account when comparing his results with ours. Indeed, in addition 
to Agc = 180°, /3qc = 0°, he assumed perfectly ecliptic orbits, i.e. with / = 0°, obtaining 
only radial and transverse components of the perturbing acceleration. Thus, his precession 
of the longitude of perihelion w redu ces to tha t of th e argument of perihelion to because 
there is no precession of the node Q. IMilgroml (120091 ) himself acknowledges that such an 
approximation is not valid for bodies like Pluto and Icarus showing high inclinations to the 
ecliptic. Actually, the quadrupolar field of X/EFE does induce a secular precession on Q as 
well, as we will see in Section [3j 

Let us, now, reason in terms of a rock-ice planetary body. By assuming for it a mass 
as large as that of Mars we have for its distance 



An Earth-sized body would be at 



?"Fie 



rpit 

'"Fie 



84 ± 9 au, 
71 ± 19 au. 



178 ±20 au, 
150 ± 40 au, 



(25) 
(26) 

(27) 
(28) 
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while a gaseous giant like Jupiter would be at 



rpit = 1.218 ± 0.135 kau, (29) 
r Fie = 1.027 ±0.274 kau. (30) 



The distance of a brown dwarf with M = 80Mj would be 



r Pit = 5.246 ± 0.583 kau, (31) 
r Fie = 4.425 ± 1.180 kau, (32) 

while an object with the mass of the Sun would be at 

rpit = 12.366 ± 1.374 kau, (33) 
r Fie = 10.430 ±2.781 kau, (34) 

Incidentally, it may be of some interest to compute the tidal parameter of Sgr A* itself to 
see if its action could be det ected from its influence on the motion of the solar system's 
planets. By assuming for it (IReid fc Brunthalerl 120041 ) M x = 4 x 1O 6 M and r x = 8.5 kpc, 
we have 

/CsgrA* = 3 x 10" 35 s" 2 , (35) 

which is 9 orders of magnitude smaller than the present-day level of accuracy in measuring 
K. 

Let us, now, abandon the direction of GC, and, consequently, the MOND scenario, and 
look at /C as a function of the ecliptic longitude and latitude of X without assuming any 
a-priori limitations on them. Let us, first, use the Pitjeva result of eq. ([!]). It turns out that 
the maximum value of /C, and, consequently, the minimum value for rx, occurs for 

A x = 18.3°, (36) 
fa = -89.9°, (37) 

i.e. perpendicularly to the ecliptic; 

/C max = (4 ± 1) x 10~ 26 s" 2 . (38) 

In this case, the distances of X, for different values of its postulated mass, are 

^Mars = 67 ± 7 au, (39) 

r Earth = 141±15au, (40) 

rjupiter = 969±107au, (41) 

fbmwn dwarf = 4. 175 ± 0.463 kau, (42) 

r Sun = 9.841 ± 1.093 kau. (43) 
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The minimum for /C, corresponding to the maximum for rx, occurs for 

A x = 182.8°, (44) 
Px = 1-7°, (45) 

i.e. basically in the ecliptic; 

/C min = (2.0 ± 0.7) x 10" 26 s" 2 . (46) 

The heliocentric distances for X are as follows 

r Ma rs = 84±9au, (47) 
r E arth = 178±19au, (48) 
Oupiter = 1.220 ±0.135 kau, (49) 

^brown dwarf 

5.261 ± 0.584 kau, (50) 
r Sun = 12.400 ± 1.377 kau. (51) 

Such results are quite similar to those obtained for GC. 



In the case of the result by iFienga et al.l (l2009bl ) of eq. (j2J), the minimum and the 



maximum of the tidal parameter occur at the same location as before. The maximum is 

/C max = (7 ± 5) x 10- 26 s" 2 (52) 



and yields 



^Mars = 57 ± 15 au, (53) 

r Earth = 120±32au, (54) 

rjupitcr = 817±218au, (55) 

'"brown dwarf 

3.521 ± 0.939 kau, (56) 

r Sun = 8.300 ± 2.789 kau. (57) 



The minimum is 



/C min = (3.5 ± 2.7) x 10~ 26 s" 2 (58) 
and the corresponding maximum distances are 

r Ma rs = 72±19au, (59) 

r E arth = 151±40au, (60) 

rjupiter = 1.030 ±0.274 kau, (61) 

^brown dwarf 

4.437 ±1.183 kau, (62) 

r Sun = 10.459 ± 2.789 kau. (63) 
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3. The action of a distant body and of MOND's External Field Effect on the 

node of a planet 

Calculating the secular precession of the node Q by means of the standard Gauss 
perturbative approach with eq. ([3]) for generic values of and of Ax, Px yields a 

non-vanishing effect. Also in this case, the exact formula is rather cumbersome; it is 

n) = — H(I, Px, Ax, n)J(I, Px, Ax, n, w), (64) 

/ x 4nv 1 - e 2 



with 
and 

in which 



7i = cos/ sin Px — sin / cos Px sin (Ax — fi), (65) 
J =(-2- 3e 2 + 5e 2 cos 2uj)Q — 5e 2 sin 2uj cos Px cos(Ax — fi) , (66) 

C? = sin/ sin /3 X + cos / cos Px sin (Ax — (67) 

Note that, in general, it is not defined for / — > 0°. If and when also the corrections AQ to 
the standard node precessions will be estimated, eq. (1641) could be used together (w) x and 
Am to constrain Ax and Px by taking their ratio and comparing it to the predicted one 
which would be independent of /C itself, being a function of Ax, Px alone. In particular, it 
would be possible to check if XqcPgc satisfy the equation 

^ = y(Xx,fc), (68) 

where 

(n) 

y(\x,Px) = -jT{± (69) 
is the theoretically predicted ratio. This wou ld also be a crucial test for the form of the 



MONDian EFE proposed by iMilgroml (120091 ). Indeed, if XqcPgc would not satisfy eq 



f l68|) . it should be rejected. The opposite case would not yet represent an unambiguous proof 
of MOND because there would still be room for the action of an ordinary planetary-sized 
body; clearly, should observational efforts aimed to detect it be infructuous, MOND would 
receive a strong support. 



4. Discussion and conclusions 

We have shown that a putative distant body X, not yet discovered, would induce 
non-vanishing secular precessions of the longitudes of the perihelion and the node of a 
known planet P of the solar system. In particular, the resulting perihelion precession would 
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be retrograde so that it would be able to explain the anomalous perihelion precession of 
Saturn recently determined from an analysis including radio-technical data from Cassini. 
An investigation of the tidal parameter of X as a function of its ecliptic longitude and 
latitude showed that its maximum value occurs for X located perpendicularly to the ecliptic, 
while its minimum occurs for X lying in the ecliptic. Accordingly, it has been possible to 
determine the present-day distance of X for different postulated values of its mass. Rock-ice 
planets as large as Mars and the Earth would be at about 80 au and 150 au, respectively, 
while a Jupiter-like gaseous giant would be at approximately 1 kau. A typical brown dwarf 
(M = 80Mj) would be at about 5 kau, while Sun-sized body would be at approximately 10 
kau. If it is difficult to believe that a main-sequence Sun-like star exists at just 10 kau from 
us, the distances obtained for terrestrial-type planets are substantially in agreement with 
theoretical predictions existing in literature about the existence of such bodies which would 
allow to explain certain features of the Edgeworth-Kuiper belt. Incidentally, let us note 
that our results rule out the possibility that the hypothesized Nemesis can be the Sun-like 
object X that may be responsible of the anomalous perihelion precessions of Saturn, also 
because, at approximately just 10 kau from us, its orbital period would amount to 1-10 
Myr, contrary to the 26 Myr periodicity in extinction rates on the Earth over the last 
250 Myr which motivated the Nemesis proposal. Moreover, our Sun-sized body X would 
not penetrate the Oort cloud which is believed to extend from 50 kau to 150 kau. The 
tidal parameter of Nemesis would be, instead, 2 — 4 orders of magnitude smaller than the 
present-day level of accuracy in measuring it (10 -26 s" 2 ). On the other hand, if our X had a 
distance of about 88 kau, as predicted for Nemesis, our result for its tidal parameter would 
imply a mass of 3OOM . 

For a particular position of X, i.e. along the direction of the Galactic Center, our results 
hold also for the recently proposed form of the External Field Effect in the framework of 
MOND in the sense that it would be able to explain the perihelion precession of Saturn 
in such a way that it mimics the existence of a body in the direction of the center of the 
Milky Way. The associated parameter q ranges from 0.2 to 0.4 — 1, while the theoretical 
predictions for various choices of the interpolating function and various values of the 
Galactic field at the Sun's location are 10" 2 < -q < 0.3. 

Anyway, further data analyses of enlarged radio-ranging datasets from Cassini by 
different teams of astronomers are required to confirm the existence of the anomalous 
perihelion precession of Saturn as a real physical effect needing explanation. 

Finally, let us note that a complementary approach to the problem consists of 
re-analyzing all the planetary data with modified dynamical models explicitly including 
also a planet X and solving for a dedicated parameter accounting for it. 
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